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J> ' Abstract. 

"^" . Various data analyses of the Cosmic Microwave Background (CMB) provide 

!TX | observational hints of statistical isotropy breaking. Some of these features can be 

^+ ■ studied within the framework of primordial vector fields in inflationary theories which 

r | , generally display some level of statistical anisotropy both in the power spectrum 

f^*i ' and in higher-order correlation functions. Motivated by these observations and the 

recent theoretical developments in the study of primordial vector fields, we develop 
the formalism necessary to extract statistical anisotropy information from the three- 
point function of the CMB temperature anisotropy. We employ a simplified vector field 
K^ , model and parametrize the bispectrum of curvature fluctuations in such a way that 

5-h ' all the information about statistical anisotropy is encoded in some parameters Xlm 

(which measure the anisotropic to the isotropic bispectrum amplitudes). For such a 
template bispectrum, we compute an optimal estimator for Xlm and the expected 
signal-to-noise ratio. We estimate that, for f^L — 30, an experiment like Planck can 
be sensitive to a ratio of the anisotropic to the isotropic amplitudes of the bispectrum 
as small as 10%. Our results are complementary to the information coming from a 
power spectrum analysis and particularly relevant for those models where statistical 
anisotropy turns out to be suppressed in the power spectrum but not negligible in the 
bispectrum. 



1. Introduction 

Cosmic Microwave Background (CMB) anisotropics have so far been analysed under 
the assumption of a statistically isotropic distribution. However, recent observations 
of "anomalous" features in the temperature WMAP maps might be interpreted as an 
indication of statistical anisotropy (although the a posteriori choice of statistics could 
make their interpretation difficult pQ). These features concern (see [2] for a recent 
review): a "cold spot" in the southern galactic hemisphere [Sill]; an alignment between 
the quadrupole and octupole modes in the temperature anisotropies [51 El [9] ; an 
asymmetry in the large scale power spectrum and in higher-order correlation functions 
between the northern and the southern ecliptic hemisphere [TUJ [TTJ [121 13] • Other 
interesting analyses look for a quadrupolar modulation of the CMB power spectrum 
[HI UHl [16], which we discuss later in more details, see Eq. (JTJ [j]. These experimental 
findings have been partially responsible for a renewed interest in models of the early 
Universe predicting some level of statistical anisotropy. 

Inflationary models where primordial vector fields play a non-negligible role are 
an interesting possibility. Models of this type have been proposed since 1989 [21] 
in an attempt to achieve an exponential expansion for the early Universe without 
having to resort to scalar fields. A similar motivation inspired more recent papers 
[221 1231 EH [251 ES]. Another possibility is provided by more general scenarios where 
vector fields do not necessarily drive inflation J2ZII2SII2SII3DIEI1I321I331E11I351I3SII3Z!. 
Some interest has also been directed to models of anisotropic dark energy, e.g. 
[331 133 EDI ED E21 S3] • See also [SI ES] for other anisotropic scenarios. 

A general prediction of inflationary models where one (or more) vector field is 
present is that the power spectrum of primordial curvature perturbations acquires a 
directional dependence on the Fourier wavevector. A widely used parametrization is the 
following [28] i 

1 + G(k) (k ■ iV X ' 

where G{k) is the (generally scale-dependent) amplitude of the statistical isotropy 
breaking induced by a preferred spatial direction N . Notice that the isotropic part of 
the expression has been factorized in the function P lS0 {k). Such a parametrization of the 
power spectrum is valid for most vector field inflationary models. A minimum- variance 
estimator for the parameter G was first built in [16] showing that an experiment like 
Planck should be able to detect a quadrupolar anisotropy in the power spectrum as small 
as 0.5% at la (see also [17]). Afterwords, in [H], ^5\ and p], the five-year WMAP 

X Let us also recall the lack of power of the temperature two-point correlation function on large angular 
scales (> 60°), asymmetries in the even vs. odd multipolcs of the CMB power spectra (parity symmetry 
breaking), both at large [T7l [18] and intermediate angular scales (see for a discussion [lj, [19] and Rcfs. 
therein). See also [20] . 

§ For a more general parametrization see [46] . 



P(k) = P lso (k) 



(1) 



temperature data were analysed in order to estimate the magnitude and orientation 
of statistical anisotropy for a power spectrum with a form as in Eq. (PQ). Assuming a 
constant amplitude G*, the authors of [IE] find G> = 0.29 ± 0.031. In [15] and later in 
[H] , the authors point out that the existence of beam asymmetries, uncorrected in the 
maps, should be accounted for and can be entirely responsible for the aforementioned 
power anisotropy (see however also the different conclusions of [16]). Anisotropic effects 
on CMB polarization and on the power spectrum of temperature fluctuations were 
also investigated in [321 EQl EU H7] and a study of a statistical anisotropic power 
spectrum for the temperature fluctuations has been presented in [52] accounting for 
anisotropies originated at the last scattering surface, in addition to the ones present 
in the primordial power spectrum. As far as constraints from Large-Scale-Structure 
analyses are concerned, Ref. [53] considered a sample of photometric Luminous Red 
Galaxies from the SDSS survey to assess the quadrupole anisotropy in the primordial 
power spectrum of the type described by Eq. (pQ). Assuming the same preferred 
direction singled out by [2], they derive a constraint on the anisotropy amplitude 
G* = 0.006 ± 0.036 (lcr), thus finding no evidence for anisotropy. Marginalizing over 
N with a uniform prior they find —0.41 < G* < 0.38 at 95% C.L. These results could 
confirm that the signal seen in CMB data is of systematic nature. However it must be 
stressed that CMB and LSS analyses probe different scales, and the amplitude of the 
anisotropy generally is scale dependent G = G{k) [28] . 

On the other hand, the study of inflationary models with vector fields has been 
also motivated by the possibility of predicting higher (compared to the standard slowly 
rolling scalar field model) levels of non-Gaussianity, combined with some level of 
statistical anisotropy in the three or higher-order correlation functions of primordial 
curvature perturbations ( "Anisotropic non-Gaussianity" ) [291 EH EU E3 EE1 EU E3 E3 
[50] (see also [61] for an overview of all of these models). 

In the spirit of [281 SHI EE] and motivated by a strong interest in anisotropic non- 
Gaussianities, we construct an estimator to extract statistical anisotropy information 
from the CMB bispectrum. We assume a primordial origin for violation of statistical 
isotropy as we investigated in [571 EH] where very general models of inflation 
with primordial vector fields have been considered including in particular possible 
contributions of non-Abelian SU(2) gauge self-interactions. For the minimal and 
simplest version of vector field inflationary models, we study a specific bispectrum 
template for which the departure from statistical isotropy is quantified by some 
coefficients Xlm- Essentially these parameters correspond to the ratio between the 
non-linearity amplitude of the anisotropic bispectrum and the isotropic one: Alm ~ 
(/nl//nl) ( see Appendix A). For such parameters we perform a Fisher matrix analysis 
to give an estimate of the potentiality of present and future experiments to put 
constraints on such anisotropic features. 

One of the motivations for our work is that an analysis of anisotropic non- 
Gaussianity can be complementary to that probing anisotropic features in the power 
spectrum, so as to offer a possible cross-check for statistical anisotropic signatures. 



There is, however, another important motivation to point out. Various models for 
isotropy violation can indeed be characterized by a negligibly small (unobservable) 
level of statistical anisotropy in the power spectrum, but they still display a relevant 
anisotropic amplitude in the bispectrum (see, e.g., [291 EH EH]). The essential reason is 
that the non-linearity parameter for the anisotropic bispectrum f^ L turns out generally 
to depend not solely on G but also on other mo del- dependent parameters which can 
compensate for a parameter \G\ <C 1. For example the case of non-Abelian primordial 
vector fields [571 EHj is particularly clear since f^ L does not even depend on G when the 
main contribution comes from the intrinsic self-interactions of the vector fields. This 
is because such self-interactions just determine the bispectrum without entering into 
the (tree-level) power spectrum. In other models f^ L ~ G 2 f(A), where /(A) is some 
functions of the vector fields involved (one can check this is the case for the models 
discussed in, e.g., Refs. [291 EIj). Therefore, even when |G| <C 1 is not observable 
through an analysis of the power spectrum, the parameters space of the models can still 
allow for a non-negligible non-linearity amplitude of the anisotropic bispectrum part. In 
fact generally it turns out that (even for unobservable level of G) | Alm| ~ |/nl V/nl 1 — 1; 
i.e the anisotropic amplitude of the bispectrum can even be as large as the isotropic part 
(see Appendix A for details). 

In this paper we will quantify what is the level required for the CMB primordial 
bispectrum statistical anisotropy to be measured, showing that there might ample room 
for a high level of statistical anisotropy in the bispectrum when \G\ <C 1. Searches for 
signatures of statistical anisotropy limited to the power spectrum would thus be blind 
to such features. 

The paper is organized as follows: in Sec. 2 we compute the three-point function 
of the temperature anisotropics harmonic coefficients, considering a simplified vector 
field model of inflation; in Sec. 3 we derive an optimal estimator for the anisotropy 
parameters; in Sec. 4 we estimate the signal-to-noise ratio; in Sec. 5 we draw our 
conclusions. In Appendix A we provide an overview of the bispectrum predictions 
of different primordial vector field models and discuss the parametrization of the 
bispectrum adopted in the paper, discussing also in some details a specific example; 
Appendix B includes a brief discussion of Bayesian analysis for complex parameters; 
in Appendix C we discuss in details the order of the Edgeworth expansion that we 
have employed; finally, Appendix D collects some useful properties of Wigner-3j and 6j 
symbols. 

2. Anisotropic CMB bispectrum 

As discussed in the introduction, several inflationary models have been proposed that 
incorporate primordial vector fields. For many of these theories, the main outcome is the 
introduction of statistical anisotropy features in the correlation functions of curvature 
perturbations £. We are interested in the three-point function, or bispectrum 

(C fycfycfy) = {2ir) 3 6^{k 1 + k 2 + k 3 )B(h, k 2 , k 3 ) . (2) 



Similarly to the power spectrum, Eq.flTJ, in these models also the bispectrum 
B(ki, k 2 , k 3 ) exhibits a dependence on the angle between a preferred spatial direction(s) 
specified by some vector (s) and the orientation in space of the wave vectors k±, k 2 and 
k%. For single vector field models a useful parametrization is (see Appendix A for a 
detailed overview of the models to which this parametrization applies) 

B(h, 4 h) = B is °(k u k 2 ) [1 + g B (k!, k 2 ) (p(h) (kx -N A y+ P (k 2 ) (k 2 ■ N A y 

+ P(h)p(k 2 ) (h ■ k 2 ) (k 2 ■ N A ) (h ■ Na) )] + 2 perms., (3) 

where N A is a preferred spatial direction, B lso is an isotropic function and p and g& are 
anisotropic coefficients which generally are a function of the moduli of wavevectors (see 
Eqs. (1A3I) through ( 1A6I) for an explicit definition of these functions). Eq.([3]) is therefore 
the sum of an isotropic bispectrum plus an anisotropic contribution. For the kind of 
models we are considering, the isotropic bispectrum ca be taken of the "local" form with 



B"'(kuh)=lhL-^- 3 , (4) 

where A is the amplitude of the power spectrum of the primordial gravitational poten- 
tial Pq>(k) ex k n ~ i taken to be scale invariant (similar definitions hold for B(k\,k^) and 
hrB(k 2 ,k 3 )). 

The bispectrum in Eq. Q will be considered as the primordial source for the bispec- 
trum of temperature fluctuations; we will build an estimator for the parameters that 
determine the anisotropic correction, treating the bispectrum amplitude in Eq. (jlj) as 
an input of our analysis (see the end of this section for a brief discussion on this) . 

A preliminary computation consists in moving from momentum to harmonic space. 

In the next few equations, we shortly review these general steps before returning to our 

specific calculation. 

The temperature anisotropies are expanded in spherical harmonics 

AT(n) ^ 

— ^ — = 2_^ a i™ Y im{n), (5) 

lm 

where the ai m coefficients are related to the primordial gravitational potential by 

a lm = Arr(-t) 1 J -^A^kM^Y^k), (6) 

A^(A;) being the radiation trasfer function. The bispectrum of temperature fluctuations 



is 



d 3 k! d 3 k 2 d 3 k 



Bt.tX = (aiimiflfe^fli.™.) = {^?{-i) h+h+h j (27r)3 (27r)3 ^A (fci)Az 2 (fe 2 )Az 3 (fc3) 

X ($(^)$(fc 2 )$(fc 3 )>y iimi (fcl)^ m2 (fc2)^ 3 (fc3). (7) 

|| For various shapes of primordial non-Gaussianity, see, e.g. |62j 
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Writing the bispectrum of curvature perturbations in the standard way as 

(C(h)C(h)C(h)) = (2n) 3 5^(k 1 + k 2 + h)B(h, 4 k 3 ) = (J) mhWhWh)) (8) 

and using 

I rf3 xe «.(fc 1+ fc 2 +fc3) = (2 7r ) 3 5 (3) (A : 1 + k 2 + k 3 ), (9) 

together with the Rayleigh expansion 

jits = ^Y,i l 3i{kx)YC m {k)Y lm {x), (10) 

lin 

we get 

BtXt 3 = (f) (~) (-i) h+h+h J dxx 2 j dhdhdhihhhfA^A^A^h) 

x E E EW^^M^M^CM [ dn k Y hm Xh)Y; im ,Cki) 

l 1 m 1 l 2 m 2 l 3 m 3 

(11) 



x J dfoYtM@)YiM@) Y te@)- 

Let us now replace our parametrization Q for the bispectrum. The anisotropic 
part of the bispectrum from ([3]) has a form of this type 

B(k u 4 fc 3 ) D £T° x g B x p^ x (% ■ %) (% • JV X ) (% ■ iV A ) , (12) 

where the indices i,j run over the three wave vectors, Pij = p(^) if i = j and 
Pjj = p(ki)p(kj) otherwise. The additional assumption we will make at this point in 
order to simplify our analysis is that the functions qb and p are actually scale-invariant. 
This hypothesis is easily met by the models we are dealing with (see e.g. [57]). 
For % = j we can use a simple expansion in spherical harmonics, e.g. 

p(fci)3B(fci,fc 2 ) (h -N A y = J2 X LMY LM (k 1 ); (13) 

LM 

the mixed (i ^ j) terms can instead be expressed in terms of bipolar spherical harmonics 
as, for instance (see, e.g., [631 IM] ) 



p(kMh)g B (k u k 2 )(h-k 2 ){k 1 -N A ){k 2 -N A )= Y, Khi Y h® Y h}LM, (14) 

h,l 2 ,L,M 

where 

{Yh^Y h } LM = y j C limil2m2 Yi imi (k 1 )Yi 2m2 (k 2 ) (15) 



mi ,7712 



and C/j^f z are the Clebsch-Gordan coefficients. 



In order to avoid any confusion, throughout the paper we will use upper-case 
indices LM to characterize the magnitude of statistical anisotropy, while the lower- 
case indices Ira will be used for the temperature anisotropies. We will be interested 
only in coefficients with L > 1, since L = = M corresponds to an isotropic monopole 
contribution. In particular, terms like those in Eq. ( TX3T) are quadrupolar contributions 
and, using the orthogonality of the spherical harmonis and the spherical addition 
theorem for the Legendre polynomials, Pi{h ■ n') = (47r/(2/ + 1)) J2m=-i Yi m (n)Y l ^ n ('h') , 
one easily finds 

*lm = ^(j>9b)Y£ m {N a )6 L 2, (16) 

15 

for L > 1. However, we will be completely general in using the expansion on the r.h.s. of 

Eqs. (ll3j) and ( fT^l) and only at the end of Section H] we will specifically apply our results 

for L = 2, a bispectrum model with a quadrupolar dependence. 

We will now write separately the isotropic contribution and the anisotropic ones 
(( 1T3|) and (FT4j) ) to B^^^. The isotropic contribution is 

B l ^2 z =Ki,i 3 G l ^X mz = Q) 0) f dxx 2 f dhdhdUhkik^A^A^A^h) 

x j^hx^ihx^ihx^ih, k 2 , ^Cn.3 , ( 17 ) 

where b^i^ is the reduced bipectrum JB5]; the Gaunt integral is given by 

G l ^t 2 m 3 = J dn & Y limi (x)Y l2 m 2 (x)Y hm3 (x) 

(2Z 1 + l)(2/ 2 + l)(2l 3 + l) ( h k h 

All I Vtl\ ?T1>2 ^3 

and B iso (k u k 2 ,k 3 ) = B iso (k u k 2 ) + B is0 (k u k 3 ) + B iso (k 2 ,k 3 ). The sums over l[ and 
m i (i = 1,2,3) were all simplified using the orthogonality condition of the spherical 
harmonics 




d^Y^MY^ik) = 5 kl ,5 mi< . (19) 

The contribution from anisotropic terms like the ones in Eq. ( TTB1) instead gives 

B hhh(A) _ ( 3 X 6 \ V^V^ A /^iWa Qhl'iL (_ 1 )h( i \h+l' 1 (_ 1 \m' 1 A 



3x6 



3 



:1 — : ) WAla/G Zi ' 2 ^ G^ , (-l) /2 (i) ?2+/ 2(-l) m 2 & z 2 



/3^6\3 ; ;l (_ 1 )fa^fa-H , a(_l)»H 6 ' 

\ 57]- / -^ ^1^ ^ ivi m 1 m 2 m3 m 3 -m 3 M v y v ' v 7 'l'a« 
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(20) 



where 

h hhh — I dxxI I dk 1 dfadfa(k 1 fafa) 2 A ll (k 1 )A l2 (fa)A l . i (fa) 

x j^fax^fax^fax) (B iso (fa, fa) + B is °(fa, fa)) , (21) 

b L } 2h = fdxx 2 f dfadfadfa{k l fafa) 2 A h {k l )A l2 {fa)A h {fa) 

x ^(M^CM^CM {B lso (h, fa) + B is °(fa, fa)) , (22) 

b hl j 3 = fdxx 2 f dk l dfadfa{k l fafa) 2 A h {k l )A l2 {fa)A h {fa) 

x j h {hx)j l2 {fax)j ll3 {fax) {B is °(fa, fa) + B is °(fa, fa)) . (23) 

Eq. (TT91) was employed for i = 2,3, whereas for % = 1 we have an integral over three 
spherical harmonics which produces the extra (compared to the isotropic case) Gaunt 
integral appearing in Eq. (BUI) . 
Finally, let us consider the anisotropic terms like those in ([1] 



jdIi h h -> I ^ 1 f ^ 1 V^ V^ \LM n LM r l'A l -A rihl'J r lA l " 

m 1 m 2 m 3 c I / / , / > H I ml m m,m~m a mi -m,m 



T T*Jl' l" ' I' l' I ll I 

LMl I m m l 1 m 1 l 2 m 2 



I .11 II \_* I f \^l I I V_^ I ff 

m i m m 1 m 2 m^ m\—m 1 m m^—m^rn 



x (_i)r»i+«4(_i)'i+'a(j)'i+'i (i)^b^ la + perms., (24) 

where 
b\il,= fdxx 2 [dk 1 dfadfa{k 1 fafa) 2 A h {k 1 )A h (fa)A h {k, 



x J A {k x x)j^{fax)j h {fax) {B is °{fa, fa) + B is °(fa, fa)) . (25) 

Notice that in the vector field models we are considering, the isotropic bispectrum 
is separable (remember Eq. (jlj)), i.e. we can introduce three functions X , Y and Z 

{fafafafB^ifa, fa, fa) = X(fa)Y(fa)Z(fa) + perms. (26) 

that give 

Khh= / dxx 2 X h (k 1 )Yi 2 (fa)Z h (fa) + perms., (27) 

b hhh = [ dxx 2 xf i (fa)Y l2 (fa)Z h (fa) + perms., (28) 

b$ 2la = J dxx 2 X l ^{fa)Y^{fa)Z h {fa) + perms., (29) 

where 

X h (x)= f dk l A h {k l )X{k l ) ]h {k l x), (30) 

Xll(fa) = j dk x ^ h {k 1 )X{k 1 )J A {k 1 x), (31) 

and so on for the other integrals in f l27|) . 



We will now further simplify our analysis by considering an anisotropic contribution 
as in Eq. f fT3|) only. We have checked that the anisotropic terms as in Eq. (I14p . for L > 0, 
contain either contributions that can be recast as the angular decomposition of Eq. ( TT3l) 
or contributions that display an angular dependence that is actually different from 
the one of Eq. ( ITS]) and which would correspond to a different anisotropic signature. 
Therefore a distinction between these two contributions can be drawn and from now on 
we focus only on an analysis of the Alms coefficients of the anisotropic bispectrum. This 
means that our results will be targeted only to the anisotropic features of the kind of 
Eq. (115]) . A similar dedicated analysis can be performed to other types of bispectrum 
anisotropic features]^ Therefore, let us write the bispectrum, given by the sum of an 
isotropic and an anisotropic part, as 

B h h l z ■■= (n, 01 en ) = B hhh ^ +R*ife*3(- A ) 

— f ( TDhhW) I J-Vi uhhh(A)LM\ \ /oo\ 

- JNL I -D mim2m3 l/jvL=l + 2_^ A «^mim 2 m 3 \f N L=l I J \ 6Z ) 

\ LM J 

where Bm^2 3 and Bmrm^rm can be read off Eqs. ( ITTj) and ( 120]) respectively. We have 
factorized the (isotropic) non-Gaussianity amplitude /jvi, defined in the traditional way 
as proportional to the ratio between the (isotropic) bispectrum and the (isotropic) power 
spectrum squared in the equilateral configuration. 

Estimators for the \lm 's will be computed in the next Section. Before proceeding 
though, it is necessary to make some further considerations. 

First of all, for the rest of the paper, we will select vector field models where 
the bispectrum of curvature fluctuations is real, which is equivalent to assuming parity 
conservation in our theories; indeed, it is easy to realize that {Q Q Q ) = (£ - £ - ( - ) , 
therefore imposing the reality of the bispectrum corresponds exactly to require that 

^kSk2^j = (C^C^C^), (33) 

i.e. 

B{-k h -k 2 , -h) = B(k u h, k 3 ) (34) 

If We have studied in some details the angular dependence of the l.h.s. of Eq. (Tl4|) . We have seen 
that a quadrupolc bispectrum anisotropy is generated, with L = 2. As an example of what we 
mentioned above, it is easy to realize that, looking at Eqs. (Tl4|) and (fT5|) . for (li = 0,h = 2), or 
for (/i = 2,^2 = 0) one just recovers an angular dependence as in Eq. (|13j) . while the other remaining 
combination {l\ = l? = 2) displays a different anisotropic signature. 

Notice also that we have found that Eq. (fT4|) (multiplied by B lso (ki,k2) appearing in Eq. (j3|)) for 
L = contains an additional isotropic bispectrum. We have carefully checked that such an isotropic 
bispectrum has a significant correlation (more than 72%) with B zso (using the correlation defined 
in [62]). In this way, when we factor out the isotropic bispectrum B tso in Eq. (j3]), the effect of 
this additional isotropic bispectrum from the second line of Eq. ((3|) is to produce with reasonable 
approximation a constant term that can be reabsorbed as described later in Eqs. (|39j) and (|40|) . 

9 



which, (using Eq. (12"U]) ). is only satisfied if 

L = even. (35) 

Also, consider the well-known property a* lm = (— l) m a;_ m , then 

l £>h h h \* _ (_-l \mi+m,2+m3 n h h h (^fil 

\m\m2m3j \ J ~mi—m2—m3' \ / 

The latter equation, together with ( |35|) . are only satisfied if 

x *lm — ( _1 ) Al-m- (37) 

Notice also that if ( |35l) and ( |37l) are true, then Bm^rn 2m3 is real. 

With Eq. ( 135]) in mind, we can now split the sum on the r.h.s. of Eq. (Tl3|) as follows 
J2 X LMY L M(k) = ^ooY 00 + Yl X LMY LM {k). (38) 

LM L>2,M 

The first term of the r.h.s. is a constant (remember that we are assuming no momentum 
dependence for the anisotropy parameters), so it can be conveniently incorporated into 
B lso through a rescaling 

B ia °^B is °(l + 2\ Q0 Y 00 y, (39) 

similarly, the anisotropy parameters would be rescaled 

Xlm -+ 1 I 9A Y ■ (40) 

With these rescalings, the L = contribution will be included in -Bmi 2 m 2 m 3 and so the 
sum in Eq. ( J32l) will be from now on intended over L > 2. 

Finally, we would like to quickly comment on the a priori determination of Jnl 
that, as previously mentioned, we will treat as an input parameter while estimating the 
Alm's. /nl can be estimated from the angle averaged bispectrum defined as 

B hl2k se £ ( k k k ) <i m3l (41) 

z — ' \ m\ rri2 ra? / 

mi,m2,m3 \ / 

to which the anisotropic part provides, as expected, a vanishing contribution 



771 1 777,2 "^3 

mi,m2,m3 * 



(42) 



(using Eq. ([20]) and I > 2). In g2J We define c Zi = V2/7+1 (i = 1,2,3) and 
n W2 ; 3 = a/(2/i + l)(2/ 2 + l)(2/ 3 + 1). This last result accounts for the possibility of 
determining the isotropic non-Gaussianity amplitude independently from the anisotropic 
one. In fact, the main reason why we perform the rescaling (139]) and fj4*0]) is exactly that 
an estimator measuring the isotropic /nl from the angle-averaged bispectrum would 
eventually also pick up the monopole contribution L = 0, and so we must take it into 
account when defining the isotropic part of the bispectrum. 

10 



3. Estimator for statistical anisotropy from the bispectrum 

In this section we will derive an estimator for the amplitudes of statistical anisotropy 
Alm- We start with the expression of the probability distribution function (PDF) for a 
general non-Gaussian case employing the Edgeworth expansion [66j [67], EHJ [69] 

(1 _V 8 d d \ e" ^9> m 9 a *4^4 C l-4^4,'5 m 5 ai 5 m 5 

1 ~ 6 ^ hmiahm2ahmi) ^,^,^ 3 + '"J (2nf^ (detCf 2 ' 

(43) 

where i — 1, 2, 3, q — 4, 5, JV p is the number of pixels of a given measurement and C is 
the covariance matrix 

Clm,l'm' = (aima} m r). (44) 

We now truncate the Edgeworth expansion at third order: one can verify that these are 
all the terms we need for the computation of the estimators and of the signal-to-noise 
ratios of the Alms coefficients, if we make the choice to work at leading order in f^L 
and also to neglect all connected correlators of order larger than three (see Appendix C 
for more details). 
Using a* m = (— l)' m a^_ m and the relations 

Ct- m /- m ' = (-l) m+m 'c fm > M , (45) 

iplm/m') = C l'rn',lmi ( 46 ) 

which also apply to the inverse matrix 

C- 1 , , = (-lr+^CT 1 , (47) 

l—m,l —m v ' I m ,lm v ' 

'C- 1 , X = CT 1 , (48) 

^ lm,l my I m ,lm' v ' 

Eq. (|43|) becomes 

P (°) = [ l + l E ( a h mi ai 2 m 2 a hm:i ) J2 ((- l T ll+m2+m3C h- m uUrn4 a h^ ( 49 ) 

li,mi=l,2,3 l q ,m q =A,5,6 

v ri— 1 „ n^ 1 n — f_1 , |™ 1 +™ 2 r _1 f~ l n 

A ^l 2 -m 2 ,l5m 5 u 'hm 5 , ^l 3 - m3i l 6m6 u 'l6m 6 \ *-j ^l 2 -m 2 ,l 3 m 3 ^ y l 1 -m 1 ,l4m4 U 'Um4 

— /_l\nii+m2fi-l C^ 1 n _ ( _-i\mi+m3/^i—l C^ 1 n I 

V -V ^/i-mi,i 3 m 3 , - / «2-m2,Uni4 U ''4m4 I x j '-'Zj-mi,^"^ Z3-m3,Um4 ^ m 4y 

g _ 2 ^i 9 ,m 9 a l4m. 4 ( ^l 4 m4,l 5 m s ai 5 m 5 
V 

(2tt) JVp/2 (detC) 1/2 
Let us now consider a diagonal variance. This choice is motivated by the fact that the 
non-diagonal part of the variance would generally introduce a correction proportional 
to the statistical anisotropy amplitude of the power spectrum (see e.g. [32]) which, as 
previously mentioned, is known to be^uite small both from experiments and in the 
theoretical models we are referring to 



"i 



+ Indeed a non-diagonal variance can be generated by anisotropics generated by various experimental 
effects (such as masking or inhomogencous noise). These effects can be treated with the linear term 
widely used in the bispectrum analyses, see, e.g., Ref. [70] . 
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Eq. 



will then reduce to 



p( n \ - [l i l ST In \ ( a *hm 1 a hm2 a hm3 izll^A A 

^{a) - j^i + - 2_^ \a hmi a hm2 a hm:i ) y — — (); 2 ; 3 d 

ii,mi=l,2,3 Zl ' 2 ' 3 ^ ' 2 



7/(1 






'^lil'i^mi—ms^ 



■l) mi 



hm2 



C h Ci 3 



~Ulil2"rni—m2 a l 3 m 3 



m2-mo, a l x m,x 

2~Yjl, m a *m C l a lrn 



{2it) Np/2 (detC) 



1/2' 



(50) 



It is easy to verify that the linear part of Eq. (j5"0|) is proportional to the monopole in 
the isotropic case. Let us for instance consider the first linear term, i.e. proportional 
to a\ m : from 5 m2 - m3 and from mi + n%2 + m^ = (which is always satisfied if the 
bispectrum is isotropic) we have m\ = 0; using the identity in Eq. (1D1I) . we get h = 0, 
so the only nonzero contribution is due to ag . A similar reasoning can be applied to 
the second and the third linear contributions in Eq. ( 1501) . This does not apply however 
to the anisotropic case, so it is necessary to retain the linear part of the expression. 



We are now ready to compute the estimator for the anisotropy parameters \lm 
using the expression for the total bispectrum, Eq. fl32|) . For small non-Gaussianity and 
for small statistical anisotropy, the PDF can be expanded as 

ln(P) = ln(l + x) - - ^2 a* lm C{ l ai m + const. 



i. m 



x 



y ~ 2 



J~]ai m Ci 1 aim + const., 



i.in 



where 

Jnl v^ 

x = > 

6 ^ 

(_l)ma 



B hhh(I) I + X^ X B hhh(LM), 

J -^mirri2m3\jNL—± ' ^_ / A iM"ui 1 m 2 m3 \jNL— i - 



LM 

(-l) mi 
'°l2h°fn2-m 3 0'l ir a 1 n n 



uiil-sVmi—mz&u 



'Iimi"l 2 m2»l 3 m3 

C h Ci 2 Ci 3 



I1W3 



W1W2 ^h^2 t>Zil>i 

Let us introduce the more compact notation 



<5/ l/2 5 mi _ m2 af ). (52) 



Urns 



a Lm, a 'lr,m n a 'Urr,., ( — 1) 2 



hh{I) I /" u 'hmi u 'i2m2 U 'i3m3 

m,2mz I/jvz,=1 



\mi 



C^Clr 



~Vhh*rn\-m: i a >l 2rn2 



C h Ci 2 Ci 3 
(-l) mi 

C h Ci 3 



C/i C/ 2 



0Z 2 /3"m2— m3 fl ;j 



oil 



( Jlll2 ( ^rn,i-m2 a l 3 m3 j ' 



d(A) __ x r>(A)LM _ U \^ f>hl2h(A)LM\ ( 

B = XlmO = -X L M 2 ,-°mim 2 m 3 l/ivz,=l ^ 



a l 1 m 1 a l2m2 a hm 3 

C h Ci 2 Ci 3 



Villi 



Cu CV 



'^l2h^rn2-mz a l imx 



\m± 



\m l 



~Olt /., mi — m^ Ci/ 



"^/l /o ^777.1 — 7T7.0 W/, 



(51) 



(53) 
(54) 



/o /-< ■ J hl3 U rni~m3 a l 2 m 2 ^ ^ °hh°m\-m2 U 'l 3mz )i 



where sum over (L, M) is understood. Notice that the operators B^ and JB^ are real 
and commute with each other. So 



ln(P) ~ / wx [5« + A LW 5^» LM ] 



/. 



NL 



(P«) 2 + 2P«A LA/ P (A)LM + {X LM B^ LM y 
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- 2 Yl a *™ C i l(llm + const - ( 55 ) 

l,m 

Let us call Xlm our estimator for the anisotropy parameters. We require Xlm to be 
unbiased and optimal estimators 

(Xlm) = Xlm, (56) 

where F is the Fisher matrix whose (LM, L'M') entry is defined as [H3I7T] (see Appendix 
B for more details) 



d 2 ln(P) \ / fd\n(P)\ (d\n(P)\ 
~^d\ LM d\l, M ,/ ~ \ \dX~LM 



Fx LM x LlM , = - ax a V = i^ T^ • (5 



/ \ d K'M'J 



It is possible to prove that, if the Fisher matrix is diagonal (which is our case as we 
show in the next section) a necessary and sufficient condition for an estimator Xlm to 
be optimal is the following (see for example [69] ) 

<91n(P) 



f Hm*lm \^lm - Xlm J ■ (59) 

We will assume that all other cosmological parameters, including J'nl, are known 
beforehand and, from there, determine the estimator for Xlm- The Fisher matrix is 

Fx LM » LM = P nl (b^^B^^ . (60) 

One can check that Eq. ( |59l) is satisfied by 

Xlm = jr^— (f NL B^ LM - f NL B^B^ M ) , (61) 

if in the expression for (d In P/8Xlm) 

^ = JnlB^ lm - f NL B^B^ LM - f 2 NL Xl, M ,B^ L ' M '*B^ LM , (62) 
oalm 

one replaces the last product with its expectation value, B^ L M B^ LM — > 
{B(a)L m jg(A)LM^ This approximation is justified in the case of a low level of sta- 
tistical anisotropy (and small non-Gaussianity). Therefore, in this limit, the estimator 
is given by Eq. (IBTjl . The first term on the left-hand side of ( IBlT) has a familiar structure 
if one considers for instance the expression for the estimator of the non-Gaussianity am- 
plitude j'nl] however, in this case, we expect the subtraction of a term proportional to 
the isotropic part (B^) of the bispectrum, since all of the information about statistical 
anisotropy is encoded in B^ LM . Notice also that the second term is not subleading 
compared to first term as the second power of Jnl might apparently suggest. The Fisher 
matrix F\ LM \ LM will be computed in the next section. 
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4. Signal-to-noise ratio 

The Fisher matrix is a very powerful piece of information since it sets a lower bound on 
experimental errors via the Cramer- Rao inequality (see, e.g., [691 E21 E2 EH E3 [76] and 
references therein). In fact our estimator ( TTB1) is optimal in the sense that the Cramer- 
Rao inequality is saturated, that is the variance of the estimator equals the one given 
by the Fisher matrix 

(F" 1 ) A . (63) 



a 



ALi\/ \ A ) \lM>* 



N, / (F - 1} 

v 'Ar M A J 



The signal-to-noise ratios for our parameters Xlm are therefore 

S \ Xlm (64) 

Let us now compute the Fisher matrix as given in Eq.( |58|) . We will compute it in the 
limit where the variance is diagonal and to leading order in /atl- Is it convenient to 
report the complete expression for the anisotropic contribution to the bispetrum from 
Sec. 3 

r( A ) - 1\ V^ uhhh{A)LM\ f a hm 1 a hm 2 a hrn 3 (-1)"% <- * 

D — g A LM / r> m 1 m2mz l/jvz,=l I n n n n n u hl3 u m2—m3 u l 1 mi 

(-l) m i # (-l) mi + \ 

~~ F< r< ^hh^m 1 -mz a l 2 rn2 TT - ^ ^hh^rn 1 -m2 a l 3 m 3 J 5 1,00) 

where a sum over (L, M) is understood. As we can see from the previous 
equation, 10 contributions arise for F\ LM \ LM : 1 from the six point function 
( a hm 1 a hm 2 a i 3 m 3 a hm i a l5m5 a l6m6 ), 6 from the four point functions {a* hmi a* hm2 a* hmz a hmA ) 
and 9 from the two point function (a* m a; 4m4 ). We provide their expressions below 

\ a hm 1 a l2rn2 a l 3 rn 3 a li , m4, a hrns a lz,Tn(il — ( Jl 1 l 4 0i 2 i 5 0i 3 i 6 rnirn4 rn2m5 m3m6 Ly ^L- \ 2 Vj \ 3 

+ ^hl2"l 3 l4."l5l6"rni-m2^m 3 m4"m5-m 6 {~ 1) ^h^h^h 

+ 13 terms, (66) 

\ a l 1 m 1 a l2m 2 a l 3 m 3 a Um 4 } = 0l 1 l 2 0l 3 l 4 rni - m , 2 rn3m4 ( — l) l Cl 1 Cl 3 

+ 0; 1 ; 3 0; 2 ; 4 m , 1 _ rn3 m2m4 ( — 1) C^C^ 
+ ^ 3 ; 2 ^ 1 ; 4 5 m3 _ m2 (5 mim4 ( — l) mj C; 1 C , i3, (67) 

( a Z* 1 mi a Um 4 ) = ^ik^mim 4 C«i- (68) 

The two, four and six point functions respectively provide the following contributions 
to the Fisher matrix 

F il P X M , = fNL {V*lm) hL'6 MM ,, (69) 

^aSw = ^ (- 18S ^m) <^"W, (70) 



and 



F £fw = &l (6SL/ + 9Ei M ) bW- (71) 
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where 

( u hhh{A)LM\ 2 
1 . ^ I -Dmim 2 m 3 I 

36 <*— ' C h Ci 2 Ci 3 

'lj'2i'3i"ll!"J2i"l3 

V 6 ! V^ (-l) m2+m5 / / 1 i 2 i 2 ( J 4)iM\ f R hhh(A)LM\ / 7 ox 

Zl,Z2,'5i "12,"»5 

The Fisher matrix is given by the sum of Eqs. (1691) . (170]) and (1711) . and it turns out 
to be diagonal. Notice that the contributions from T, b LM cancel out. The variance of the 
estimator becomes 

-^ = I(3Q + 6.M) LM , ( 74 ) 



a\ 6 



where 



Q V (gi + 1) (2/ 2 + 1) (2/ 3 + 1) (2/; + 1) . 2 ( h+ / i ) . 

{4x)*C h C h C k % \ 



hhl'jl-i 



h 



h h h \ ( k h L 
/ I 



'2 



KA v- (2/i + 1) (2k + 1) (2/ 3 + 1) (2/; + 1) (2l 2 + 1) .( h+l ' 1+h+l ») V 




'2 k \ I k k L \ I I2 ( 2 





M / V 

which are independent from the multipole index M. The expression of the bispectra 

(b[\ hl3 ), (b\\ hh ), (6,^, s ) are given in Eqs. 02123}. Equations (EQ) through (US} were 
employed to derive Q and M.. 

We will now evaluate ( 175]) and (1TB]) for L = 2. The triangle inequalities involving Z x 
and / x are such that l x G [k — 2, l\ + 2]; similarly / 2 £ [^2 — 2, / 2 + 2], therefore we expect 
the approximations l x ~ l\ and l 2 ~ Z 2 to be valid in Eqs. ( 175]) and ( 176]) for / l7 Z 2 3> 1- 

In order to evaluate the sums over lx, Z 2 and Is, we will assume l™ ax ~ l™™ 10 — 
jmax ^j ^ w ] 1 i c ] 1 j s the maximum multipole to be probed by a given experiment. Under 
these approximations the various bispectra, Eqs. ( 12TH23]) . turn out to be all equal and 
related to the reduced bispectrum bin given in Eq. (TTTj) . Using Eqs. ( ID 6 1) to ( 1D8I) listed 
in Appendix D, we derive 

Q^^M-Ahufi 1 1 : \ [ : : : 1 * ( — :;— i ^^, w 







'2 x 0.36 x 0.125 



(78 
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where we replaced the sum over li, l 2 and 1% by a factor I 3 . 
The final result for the variance of the estimator is 

1 _ __ l 5 h 



0-05 - "'" . (79) 



Using the numerical estimates of l 2 Ci ~ 6 x 10" 10 and b m ~ (2/3) x 2.4 x 10~ 17 f NL l~ 4 
from [65], we finally get 

O-i/iVL (t^) ■ (80) 



a A2M — V2000 

Notice that in deriving Eqs. (ITT]) and (ITS!) under the approximations explained above, 
the sign of the contribution Q is essentially determined by the sum over l t of i 2( - ll+l i\ 
Taking l x ~ li ~ l 2 — l 3 ~ Z, all the other contributions are factored out of this sum 
and the net result is just a +1 contribution as in Eq. ( ITTj) . A similar argument holds for 
the sum over l x and l 2 in Eq. ( ITS"]) . By performing numerically the sum over l { (keeping 
all the bispectra equal to bm) and then summing over h,l 2 ,h, we have checked that 
the sign and the order of magnitude of our results are correct. In fact, in this way, we 
are able to give an improved computation of the signal-to-noise ratio which does not 
underestimate the contribution from the sums over l t . We find 

QA f"L (t^) ■ (81) 



a X2M J "^V2000, 

For a (full-sky) experiment cosmic variance limited up to l max = 2000 the secondary 
effects from lensing are expected to be still subdominant in the C/'s appearing in 
the denominator of ( 1T5]) and ( ITBl) [651 ITTj . In this case, if we take the central value 
for the local non-Gaussianity amplitude Jnl = 32 (the constraints from [T8] indicate 
— 10 < Jnl < 84, 95% C.L.), then we see that at la the variance for the amplitudes of a 
quadrupolar statistical anisotropy in the bispectrum A2A/ turns out to be <J\ 2M — 0.08. 
If we consider l max = 1500 we find cr\ 2M ~ 0.1. This value is roughly representative of an 
experiment like Planck which is signal-dominated up to l ma x — 2000. This means that 
an experiment like Planck could be sensitive to an anisotropic bispectrum amplitude of 
this type as low as 10%. 

5. Conclusions 

In this paper we built up an optimal estimator for statistical anisotropy parameters 
from the CMB bispectrum. We considered models where statistical isotropy is broken 
in the early Universe by the existence of a preferred spatial direction (TV) due, for 
instance, to a primordial vector field. The vector field fluctuations can affect, by various 
mechanisms, the correlation function of curvature fluctuations, by introducing some 
degree of statistical anisotropy. We parametrized the expression of the anisotropic 
primordial bispectrum as Eq. ([3]) and derived, from there, the bispectrum of the 
coefficients of the temperature anisotropics harmonic expansion, Eqs. (fTTj) . (120]) and 
(124]) . In particular, we focused on the simplest possible parametrization of the primordial 
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three-point function, by considering terms of the bispectrum that are proportional to 
(k ■ N) 2 . We also assumed ^-independent anisotropy parameters. For these parameters, 
we computed optimal estimators, Eq. ( 16T1) . and their signal-to- noise ratios finding, for 
a quadrupolar anisotropy, a sensitivity for a Planck-like experiment to an amplitude 
of about 10%. We have pointed out that such an analysis is particularly relevant for 
those models that predict a negligibly small amplitude of the anisotropy in the power 
spectrum and, at the same time, a potentially high level of statistical anisotropy in the 
bispectrum. In Appendix A we have given some general arguments about this point, 
and we have considered as an example some details of the vector curvaton model. But, 
this can be quite a generic situation (as one can check looking at the results of, e.g., 
[29| |3T| [57]; see also the recent work [79] on the anisotropic CMB bispectrum in the 
vector model proposed in [29]). For these models analysis focusing only on the power 
spectrum could result completely blind to such bispectrum anisotropic signatures. 

Our analysis is valid for a (full sky) cosmic variance limited experiment up to 
Z max ~ 2000 (without accounting for polarization). We have seen that an experiment 
like Planck could be sensitive to a quadrupolar anisotropic bispectrum amplitude as low 
as 10%. Looking forward, we can envisage that the polarization information can help 
in improving the sensitivity in a similar way as it has been realized for the search of the 
/nl non-Gaussianity [TTJ EHl [HH E21 [83]. Of course a realistic search will need to account 
for various systematic and foregrounds effects (e.g. see Refs. [H LT31 LT5], [T6l SHI E3]) - Also 
lensing of CMB needs to be taken into account in order to avoid possible anisotropic 
effects (see, e.g. [EH] and Refs. therein). 

Our work could be extended in several directions. For instance, it would be 
interesting to study the anisotropic corrections expressed in terms of bipolar spherical 
harmonics (Eq. f[T4l) ) and their products, so as to generalize our calculation to all 
known vector field models (thus also including non-Abelian models). Also, a correlation 
between the bispectrum and the power spectrum statistical anisotropy as well as a 
k-dependence for the anisotropy parameters could be investigated. 
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Appendix A. Primordial vector field models of inflation and the 
bispectrum. 

The bispectrum of curvature fluctuations has been computed in different inflationary 
models with primordial vector fields; in this paper we specifically refer to the results 
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from [2H] (a hybrid inflation model with an Abelian vector field coupled to the waterfall 
field), from [56] (about a single vector field both in curvaton and hybrid inflation 
models) and from [57] (on multiple vectors, in particular forming an SU(2) triplet, 
both within the curvaton model and in standard inflation). A compact expression that 
well summarizes all of the "Abelian" results (in particular characterized by only one 
direction of anisotropy N) is 



B(k 1 ,k 2 ,k 3 )=B* so (k 1 ,k 2 



rlkt-N 

n ( k\ ■ k 2 



(k 2 -N) +e(k 1 -N) (k 2 - N^ 



k 1 -N)[k 2 -N 



2 perms., 



(Al) 



where B tso is a function of the wave vectors lengths (/c; = \ki\, i = 1,2,3), that also 
depends on the parameters of the specific model; the same applies to the functions T, 
A, 9 and Q. 
Similarly, the power spectrum has the following form 

f G{k) (k-N) . (A2) 



p{k) = P lso (k) 



As mentioned in the introduction, various data analysis have been performed for a power 
spectrum as in (1A2j) to find out what the orientation of the preferred direction should 
be and to set some upper bounds on the amplitude G. No such analysis has yet been 
performed for the bispectrum. 

Considering a non-Abelian 577(2) model we can write an even more general 
parametrization of the bispectrum (taking for simplicity three vectors of the same length 
N A = \N a \ (a = 1,2,3)) 



l + # B (fci,A; 2 )(p(fci)5](fci-^ c 

a 

+ p(ki)p(k 2 )ki-k 2 Y,(h-N c 



+ p(k 2 )J2(k2-N a ] 

a 

k 2 ■ N a ) )1 +2 perms. 



B(k 1 ,k 2 ,k 3 )=B lso (k 1 ,k 2 



+ g 2 c Hl J2 ^i' fc 2, k 3 )I n {k ■ kj, N a ■ N b , k t ■ N a ) (A3) 

where the SU(2) index a runs from 1 to 3; the i,j indices in the argument of the I n 
functions label the three (unit) wave vectors. The isotropic part of the bispectrum is 
determined by the contributions of both the inflaton and the vector fields, i.e. 



B is °(k!, k 2 ) = NlN^P^k^P^h) + NlNAAP+^P+ih 



(A4) 



where the power spectra of the fields are weighted by the derivatives of the e-foldings 
number with respect to the different fields, e.g. N a = (dN)/(dA a ). Typically N a oc A a . 
Also we set N A A^ i:i = {d 2 N)/(dA a i dA b j ). Here P+ = (1/2)(P R + P L ) is the mean power 
spectrum of the two transverse polarization states of the vector field. Notice that the 
structure of the first two lines include the Abelian case (1A1|) . whereas the last line of 
( 1A3I) is specific to the non-Abelian case since it includes the contribution from SU(2) 
interactions (g c is the £77(2) coupling constant). 
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We report the explicit expressions for the functions appearing in (1A3j) 

„(*,, fe) = : | 1 + ( - ) | ^ j p +(fci) p +(fe) I < A5 ) 




**> B W ■ (A6) 

where Pj(A;) is the power spectrum of the longitudinal component. In the same context, 
the power spectrum is 



P(k) = P iso (k) 



1 + G(k) ^(k-N* 



(A7) 



where 




G = p(k)gp(k)=p{k)ll+[ T Z\ -M . (A8) 



P iso (k) = N^P^k) [1 + j3{P + /P^)\ and we introduced the parameter (3 = (N A /N^f 
(the reader can refer to [57] for more details about the quantities introduced in these 
equations). Notice the similar structure for the power spectrum and bispectrum related 
quantities, e.g. Eqs. flSTJ-O and Eqs. (fA8l)-( fA5|) . 

For most vector field models, it is correct to set p ~ 0{\) and P^ ~ P + . Then 
generally one finds that 

p (vector) 

G ~ g P ~ -*- (3 , (A9) 

where P^ ve ' - or - 1 is the (isotropic) contribution of the vector field(s) perturbations to the 
total (isotropic) curvature perturbation P^ and the last step is valid for <C 1. Similarly 
for the bispectrum 

p (vector) 

9B ~ ^tj^ , (A10) 

5 C 
corresponds to the ratio between the isotropic contribution from the vector field to the 
total isotropic bispectrum (see Refs. [57J [58] for further details). No data analysis has 
been performed for the power spectrum if more than one anisotropy direction is involved; 
however, we will safely assume that the existing bounds are also valid for a multi-vector 
case, so G <C 1 and /3<1. 

In the language of the angular decomposition (fT3|) . the coefficients Alm turns out 
to be 



,(A) ,(A) 

JNL _ J] 

Al) Al)<j> , 



A LM ~ 9bF(A) ~ ^ = ,,u ' .aw^ > ( Al1 ) 



?(!) j. (1)0 , j.(I)vector 



so that they represent the ratio between the amplitude of the anisotropic (A) bispectrum 
from the vector fields to the total isotropic (I) bispectrum (which generally receives 
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contribution from both the inflaton and the vector fields). Here F(A) account for 
some angular configurations which is model-dependent. Just because of this angular 
configuration dependence, generally it turns out that \f^ L | < | /nl C ° r | , so that one 
expects |Alm| < 1- 

5.1. A worked example 

As an example, consider those models where the contribution to the bispectrum from the 
inflaton field is negligible (as for the simplest scenarios). In this case /^l m ^ ( 1- dAJJJ) 
is negligible and it is easy to obtain from Eq. ( 1A4I) 

,(/)vcctor _ ,(/) ^ N\N AA _ 1 N AA q2 N AA ( \-[0\ 

iNL — JNL — ( M 2 , t\t2\2 ~ n , Q\2 1\t2 ~ P t\t2 " \ ALZ ) 



W + N\f (1 + ft* N\ " AT 



A 

In the case of the vector curvaton model, for example, one finds that (N AA /N A ) ~ fi^ 1 , 
with Q A the density parameter of the vector curvaton at the epoch of the curvaton 
decay, so that 

/£ = £■ (A13) 

According to Eq. (jAlip the anisotropic non-Gaussianity amplitude /^ L differs from the 
isotropic one by the angular configuration dependence encoded in F(A), which generally 
is |F(A)| < 1 (this argument agrees with the explicit formulae Eqs. (34) and (35) of 

0U). 

At this point we can answer two important questions. First, given that the 
amplitude of the (isotropic and anisotropic) non-Gaussianity depends on G 2 ~ /? 2 < 1, 
is /jvjl suppressed? This is relevant given our result Eq. (ISIj) . Second, does a small 
(non measurable) anisotropy amplitude (G) of the power spectrum imply a small (non 
measurable) amplitude of the anisotropic bispectrum (Alm)? The answer to both 
question is actually negative. For example if one takes G ~ 10 -2 then, in order to 
have /nl ~ 30, we require Q A ~ 10~ 5 — 10~ 6 , which well satisfies the requirement that 
in these models the vector field must remain subdominant (Jl A <^C 1). Notice that this 
also satisfies the condition Vt\ > f3P^ ~ 10~ 10 f3 which comes from the requirement that 
|(L4/A| < 1 in order to avoid a non-Gaussian perturbation to dominate (see [31]). [j 
This answers also the second question, since the anisotropic non-Gaussianity amplitude 
/n L differs from the isotropic one by the angular configuration dependence |F(A)| < 1. 
Therefore it is certainly possible that /^ L ~ /nl an( ^ therefore to have |Alm| ~ 1- Our 



result ( 18 ip on the other hand shows that in principle an experiment like Planck can be 
sensitive to much lower values of Alm- 

The discussion following Eq. ( 1A8[) has been developed with the Abelian 
contributions to Eq. (1A3[) in mind. For the non- Abelian contributions the amplitude of 
the bispectrum is totally independent w.r.t. the power spectrum anisotropy amplitude 

* Recall that the vector curvaton contribution to the curvature perturbation is C,a — 0, A (5A/A) and 
that Pl so = P c (l+/3). 
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and the parameter space of the theory is such that a comparable amplitude is possible 
for fx L and f N [ (see [57] for a direct comparison between these contributions) . From the 
expressions of the anisotropy coefficients I n appearing in the non-Abelian contribution 
[57] it is easy to check that, in order for these contributions to be accounted for in the 
computation of the CMB bispectrum, one would have to resort to products of bipolar 
spherical harmonics. 

However, in order to keep our calculations as simple and straighforward as possible 
in this first paper on statistical anisotropy estimators from the bispectrum, we focused 
on the Abelian case, considering the Abelian version of the first two lines of (1A3J) (given 
in Eq. ©). 



Appendix B. Fisher matrix/Bayesian analysis for complex parameters. 

Let us derive the Cramer-Rao inequality for complex parameters following the 
computation for the real case reviewed, for instance, in Sec. 3.1 of [7_I]. Given two 
parameters Ai and A2 and their unbiased estimated values X\ (i = 1,2), the following 
inequality 

(Bl) 



( (AAi + /3AA 2 ) (AAi + /3AA 2 )* ) > 0, 

(where AAj = A, — \) must hold for any value of /3. In particular, for = 
-(AAiAA2)/(|AA 2 | 2 ), Eq. (jHU) provides the Schwarz inequality 

(|AA 1 | 2 )(|AA 2 | 2 )>|(AA 1 AA;)| 2 , (B2) 

which we can use in the derivation of the Cramer-Rao inequality as follows. Consider 
for simplicity one parameter only A. By definition 

(AA) = f A\P(x,X)dx = 0, (B3) 

for some data x. Let us derive the previous equation w.r.t. A 

/ (AA) (^jpj dx + J Pdx = 0, 

then 
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where we used dP/dX = P (dlnP/dX) and J Pdx = 1. From the Schwarz inequality 
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we get 

d\nP\ fd\i\P 



/<9 2 lnP 



OX J V d\* J I \ \dXdX 
which, after replacement in Eq. (IB 6 1) , gives 



IAAI 



> 



1 



d 2 lnP 
dXdX* 



where F is the Fischer matrix for A. 



Appendix C. Higher order terms in the Edgeworth expansion. 



The Edgeworth expansion is given by 
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(C3) 
(C4) 
(C5) 

(C6) 



and so on for the rest of the series (the subscripts c in the previous equations label the 
fully connected correlators). 

We want to perform our computation of the estimator and of the Fisher matrix for the 
A^ms coefficients to leading order in powers of /jvl! this implies that we can immediately 
drop all the terms in the Edgeworth expansion with coefficients that are proportional to 
(fNL,) k , k > 3 (see e.g. Eqs. (!60l) - (l6T|) ). Moreover, we choose to neglect the four-point 
function and all higher order (connected) correlators. With these premises, the only 
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terms left in the Edgeworth expansion are the first line of Eq. (jCip and the second line 
of Eq. ([C5D 
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The next step is to compute the derivatives of the Gaussian PDF; for a diagonal variance 
we find 
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For small non-Gaussianity, Eq. (15TT) then reads 
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was also defined in Eq. ( 152]) ). Let us derive Eq. (ICIOP w.r.t. A 
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The derivative of the first bracket on the r.h.s of the previous equation is given by the 
r.h.s. of Eq. ( 16211 ; the derivative of the second bracket is 
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We can easily show that the y term in Eq. (1C10|) does not affect the estimators of the 
Xlm coefficients: Eq. ( 162]) becomes 



|^ = /vl5 (A)LM - f NL B^B^ LM - f NL Xl, w B^ L ' M '*B^ LM + -^-(y); (C14) 

the contribution from the derivative of y is null once one replaces it with its expectation 
value, exactly like it was done for instance for the term B^ L M B^ A > LM (see comment 
after Eq. (162]) ). In fact, it is easy to check that (Sji 1 *^^^) = at zeroth order 
in /tvl and neglecting the connected four and six-point functions. 
The contribution of the y term to the Fisher matrix is 

d 2 y 



Fa "' w d " \dxt, M ,dx LM ) (C15) 

- IML V 17 -\\M' uhhh(LM)\ r>kkML'-M')\ ] ( T^hhhkhk \ 

oc / j V / raim2m3 I j 'n l — 1 rrnmsme 1/jVL— 1 \m1m2m3m4m5m6 / i 

kmi 

therefore the Fisher matrix has no leading order contributions other than the ones al- 
ready derived in Sec. (3). 

We have thus shown that, if one neglects all connected correlators beyond third 
order and wants to keep the leading order contributions in f^i to the Fisher matrix and 
to the estimator, the Edgeworth expansion can be truncated as in Sec. (3). 

Appendix D. Useful relations involving the Wigner 3-J symbols. 

We will list the identities involving the Wigner 3-J symbol that we used in our 
calculations 
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